AD  A  1  3  9 1 4  4 


NUSC  Technical  Report  6885 
11  January  1984 


Extension  of  Finitely  Conducting 
Earth-Image-Theory  Results 
To  Any  Range 

Peter  R.  Bannister 

Submarine  Electromagnetic  Systems  Department 


Naval  Underwater  Systems  Center  _.I/W 

Newport,  Rhode  Island  /  New  London,  Connecticut  Ol 

,ELECTE' 

'MAR  1  9 1984 


Approved  (or  public  release;  distribution  unlimited. 


84  03  16  067 


•mo  pm  c  pnpv 


Preface 


This  report  was  prepared  under  NUSG  Project  No.  A59007,  “ELF 
Propagation  RDT&E”  (U),  Principal  Investigator,  P.  R.  Bannister  (Cede  3411), 
Navy  Program  Element  >io.  1 1401 N  and  Project  No..  X0792-SB,  Naval  Electronic 
Systems  Command  Communications  Systems  Project  Office,  D.  Dyson  (Code  TME 
1 10),  Program  Manager  ELF  Communications,  Dr.  B.  Kruger  (Code  PME  1 10-X1L 

The  analysis  and  write  up  of  this  report  was  performed  while  the  author  was 
occupying  the  Research  Chair  in  Applied  Physics  at  the  Naval  Postgraduate  School, 
Monterey,  CA.  The  author  would  especially  like  to  thank  Professors  Otto  Heinz 
and  John  Dyer  and  Dean  Bill  Tolies  for  recommending  him  to  occup/  >.n:s  post  and 
NAVSEA  tCode  63R)  for  sponsoring  the  Chair- 

The  Technical  Reviewer  for  this  report  was  Anthony  Bruno 


Reviewed  and  Approved:  1 1  January  1984 


D.  F.  Dence 

Head,  Submarine  Electromagnetic 
Systems  Department 


The  author  o!  this  teport  is  located  ai  ihe 
New  London  Laboratory.  Nava!  Underwater  Systems  Center, 
New  1  ondon,  Connecticut  06320. 


REPORT  DOCUMENTATION  PAGE 


RE*D  I'STRICTIONS 
BEFORE  CORPLETIM;  FORM 


1  WMMTNUMKR 

TR  6885 

4.  TITLE  fmkntUt 


GOVT  ACCtt&lON  NO 

\q..A/3‘7(\ 


EXTENSION  OF  FINITELY  CONDUCTING  EARTH- IMAGE- 
THEORY  RESULTS  TO  ANY  RANGE 


7  AUThOAmi 

Peter  R.  Bannister 


I.  7WOM «MQ  OMOAMZATION  NAUC  A  NO  AOOM3S 

Naval  Underwater  Systems  Center 

New  London  Laboratory 

New  London,  C'l'  06320 _ 

it  coNTMOuma  omcc  mamc  ano  aoomss 


1*.  MOWTOMNO  AGENCY  NAMC  4  AOCMSS  Uf  difftrm  frmm  Cam trmUimf  Off lew 


1  AfCmCNT't  CATALOG  NOMMA 


».  type  or  Nnm  4  pcmoo  covcmo 


4.  mwiww  owo.  wow  www 


4.  CONTAACT  OA  OAAWT  NUMMAlaf 


10.  OAOOAAM  ELEMENT.  MOJCCT.  TASK 
AAEA  4  WOKK  UMT  MUMEStt 


12.  MFOHTOATi 

11  January  1984 

12.  MMR  OK  PAOES 


IS.  StCUMTV  CLASS.  (Ail 

UNCLASSIFIED _ 

lie  OCCLASSWCATKMIOOWMOHADMa 

schedule 


14.  OtSTAWUTiOM  STATEMENT  taf  cAu  «>pMl 


Approved  for  public  release;  distribution  unlimited. 


■  OUT*  SUTTON  STATEMENT  m/  tfc.  , 


14  SUPPUEMENTAKY  NOTTS 


*  mtwrw  us  Stock  JO.  tf  different  frmm  Krfmftt 


It.  KCY  WONOS  (Cmmmmm  mm  rmmmmm  -Am  if* 


fmmd  idammtfy  by  btmek  1 


Air-to-Air  Propagation 
Electromagnetic  Fields 
Horizontal  Electric  Dipole 


Horizontal  Magnetic  Dipole 
Vertical  Electric  Dipole 
Vertical  Magnetic  Dipole 


20.  ABSTRACT  iCmmnmwm  mm  rtwfM  aid*  tf  manammry  mmd  •drtmfy  by  kimrk  mmmtbw 

"  Finitely  conducting  earth- image- theory  techniques  have  been  employed 
to  determine  new  formulas  for  the  electric  and  magnetic  fields  produced  by 
the  four  elementary  dipole  antennas  for  the  air-to-air,  surface-to-air,  air- 
to-surface,  and  surfacc-to-surfacc  propagation  cases.  The  only  restriction 
on^thc  use  of  these  formulas  is  that  the  index  of  refraction  be  large  (i.e., 
Jn  |  10).  They  are  valid  at  any  frequency  and  at  any  range  for  the  flat- 

earth  case.  These  formulas  reduce  to  previously  derived  results  when  either _ 


,^1473 


■tl  j. 


~*a j/ije.  lo-j-uj. 


,r{  A'  t  J,  =.  I  0  ^ 


20.  (Cont'c!) 


(1)  the  measurement  distance  is  much  less  than  a  free-space  wavelength,  (2)  the 
Sommerfeld  numerical  distance  is  small,  or  (3)  the  measurement  distance  is  much 
greater  than  an  earth-skin  depth. 

In  terms  of  computer  time,  these  new  formulas  can  be  evaluated  in 
fractions  of  a  minute  compared  with  hours  for  the  complete  numerical  eval¬ 
uation  of  the  exact  Sommerfeld  integrals. 

These  formulas  are  intended  to  supplement  the  author's  recently  derived 
subsurface-to-subsurface  and  air-to-air  propagation  formulas. 
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EXTENSION  OF  FINITELY  CONDUCTING  EARTH- IMAGE 
THEORY  RESULTS  TO  ANY  RANGE 

INTRODUCTION 


During  the  past  several  years,  finitely  conducting  earth-image  theory 
techniques  have  proved  quite  useful  in  determining  the  quasi-static  fields  of 
antennas  located  near  the  earth's  surface  for  both  single-layered  and  multi¬ 
layered  earths.  (For  detail  references,  see  Bannister . 1 »2)  The  quasi-static 
range  is  defined  as  that  range  where  the  measurement  distance  is  much  less 
than  a  free-space  wavelength. 

Physically,  the  essence  of  the  quasi- static-range  finitely  conducting 
earth-image  theory  technique  is  to  replace  the  finitely  conducting  earth  by  a 
perfectly  conducting  earth  located  at  the  (complex)  depth  d/2,  where  d  =  2/yj 
and  Yj  =  [iujugtaj  +  iuej)]1/2  is  the  propagation  constant  in  the  earth  (see 
figure  1  for  the  image-theory  geometry).  Analytically,  this  corresponds  to 
replacing  the  algebraic  "reflection  coefficient,"  (uj  -  X)/(ux  +  A),  in  the 
exact  integral  expressions  with  exp(-Ad) where  X  is  the  variable  of  integra¬ 
tion.3  For  antennas  located  a*  or  above  the  earth's  surface,  the  general 
image-theory  approximation  is  valid  throughout  the  quasi-static  range. 1>2 

Recently2*4  we  have  shown,  for  horizontally  polarized  sources,  that 
finitely  conducting  earth-image  theory  techniques  are  not  limited  to  the 


IMAGE 


Figure  1.  Image-Theory  Geometry 
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quasi-static  range  alone.  That  is,  by  replacing  the  horizontally  polarized 
algebraic  "reflection  coefficient,"  (Uj  -  u0)/(Uj  +  Uq),  with  exp(-u0d),  we 
demonstrated  that  finitely  conducting  earth-image  theory  techniques  can  be 
utilized  at  any  range  from  the  source.  MohsenS  has  validated  and  extended 
these  results  to  include  higher-order  terms  that  correspond  to  multiple  images 
at  the  same  location.,  Mahmoud  and  Metwally,6  employing  discrete  and  discrete- 
plus-continuous  images,  have  computed  satisfactorily  the  change  in  the  input 
impedance  of  a  vertical  magnetic  dipole  (VMD)  due  to  the  presence  of  the  earth 

We  have  also  recently  shown'.**  that,  for  small  Somraerfeld  numerical  dis¬ 
tances,  nearfield  and  farfield  range  finitely  conducting  earth- image  theory 
techniques  can  also  be  employed  for  determining  the  fields  produced  by  hori¬ 
zontal  electric  dipole  (HED)  and  horizontal  magnetic  dipole  (HMD)  antennas 
(which  are  a  combination  of  vertically  and  horizontally  polarized  sources. 

It  is  the  purpose  of  this  report  to  extend  the  use  of  finitely  conducting 
earth-image  theory  techniques  to  any  range  and  to  present  new  formulas  for  the 
electric  and  magnetic  fields  produced  by  the  four  elementary  dipole  antennas 
for  the  air-to-air,  surface-to-air,  air-to-surface,  and  surface-to-surface 
propagation  cases.  The  only  restriction  on  the  use  of  these  formulas  is  that 
ii2]  >_  10,  where  n  =  Yj/Yo-:  They  are  valid  at  any  frequency  and  at  any  range 
for  the  flat-earth  case..  These  formulas  reduce  to  the  author's  previously 
derived  results  when  either  (1)  the  Sommerfeld  numerical  distance  is  small, 7* 8 

(2)  the  measurement  distance  is  much  less  than  a  free-space  wavelength , 1 >2  or 

(3)  the  measurement  distance  is  much  greater  than  an  earth-skin  depth. 9 

In  this  report,  the  four  elementary  dipole  antennas  [vertical  electric 
dipole  (VED) ,  VMD,  HED,  and  HMD]  are  situated  at  height  h  (h  0)  with  respect 
to  a  cylindrical  coordinate  system  (p,$,z)  and  are  assumed  to  carry  a  constant 
current,  I.  The  axes  of  the  VED  and  HED  (of  dipole  moment  p)  are  oriented  in 
the  z  and  x  directions,  respectively,  while  the  axes  of  the  VMD  and  HMD  (of 
dipole  moment  m)  are  oriented  in  the  z  and  y  directions,  respectively.  The 
earth,  which  is  assumed  to  be  a  homogeneous  medium  with  conductivity  Oj  and 
dielectric  constant  Cj  (=  ere0),  occupies  the  lower  half-space  (z  <  0)  and  the 
air  occupies  the  upper  half -space  (z  ■>  0) .  The  magnetic  permeability  of  the 
earth  is  assumed  to  equal  v.Q,  the  permeability  of  free  space.:  Meter-kilogram- 
second  (MRS)  units  are  employed  and  a  suppressed  time  factor  of  exp(iut)  is 
assumed. 


AIR-TO-AIR  PROPAGATION  DERIVATION  PROCEDURE 


As  an  example  of  our  derivation  procedure,  consider  an  HED  source.  When 
h  and  z  are  >0,  the  Sommerfeld  integral  expressions  for  the  HED  Hertz  vector 
are^-O-l^ 


nx  = 


p 

e'Y°R0 

4iti(jje0 

R0 

>'yori  rVuo<=+h) 


♦  2 


f  -Mz+i 

/  - - J0(Xp)XdX 

X  U1  +  uo 


(1) 


and 


2 


TR  '.,85 


=  £ 


cos  J  d  f  2tui  ‘  u0)  -uc(z+h) 

- x  —  /  -5 - 5 — -  e  0  J_(Ap)AdA 

7Ilueo  *>I  >tu0  +  *0U1 


where 


4Ttiue 


P2  ♦  (z  -  h)2, 
>2 


(2) 


R-  =  p^  +  (2  +  h)  , 


.  xz  +  n> 

u2  =  A2  +  Y2, 

y20  =  _u)2uoeO  ’ 

,2 


and 

Yj  =  iu.u0  Co  x  +  iuCj)- 


From  equations  (1)  and  (2),  utili2ing  the  identity  (u,  -  uQ) (Uj  +  uQ)  = 
Y2  -  Yq.  we  have 
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For  in2 |  >>  1,  we  have  shown  that2  4  >7,8 
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where 


d  -  2/yj  . 

If  we  use  equations  (1)  and  (5)  and  Soraraerfeld’s  integral, 
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where  R(,  =  p-'  +  (d  +  z  +  h)2.  This  equation  is  valid  at  any  range  from  the 
source. 

Since  yJ/Yi  =  A2  =  1/n2 ,  equation  f 3 J  can  be  rewritten  as 


where 


:t  P  cos  4>  „ 
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Since  |n2|  >>  1  (jA2|  <<  1),  we  can  set  the  function  u2  in  the  second 
term  of  equation  (1A)  equal  to  Yj>  the  propagation  constant  in  the  earth, 
lheiefore, 


A2Uj  -  A2y j  =  y0a 


(13) 


and 
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(14) 


Since  JyqAI  <<;  1>  the  integral  P  will 
>>  1.  Waitl0,13  has  shown  that,  when  jn2j 


be  of  importance  only  when  | Yq^i I 
»  1  and  IyqRJ  »  1, 
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F(w)] 


(15) 


where 

F(w)  -  1  -  i(irw)  1/2e“werf c(iw-1/2)  (16) 

is  the  Sommerfeld  surface-wave  attenuation  function,  sin  1J1.  =  (z  +  h)/R1,  and 
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Y  R 
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2 — (sin  4,'j  +  A)‘ 


(17) 


is  the  Sommerfeld  numerical  distance.  For  small  numerical  distances  F(w)  -  1, 
while  for  large  numerical  distances  and  negative  arguments,  F(w)  -  -l/(2w). 

For  |n2j  >>  1,  the  Fresnel  reflection  coefficient  for  vertical  polari¬ 
zation  reduces  to 


sin  j  -  A 
sin  +  A  ' 


Since 


i  -  r, 


ii 


2  sin  iji  i  +  A 
equation  (15)  can  be  rewritten  as 
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Therefore,  from  equations  (14)  and  (20), 
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where 


[1  -  F  (w)  ] 


(^•(^) 


(22) 


F(w)  = 


sin  +  AF(w) 
sin  sj/j  +  A 


Thus,  from  equations  (9)  and  (21), 
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(23) 


Another  factor  that  we  will  encounter  in  the  derivation  of  the  field- 
strength  components  is  the  factor  B,  which  is 


E  =  sin  <ji,  -  AA 

sin2  ip1  -  A^F(w) 
sin  ij/1  +  A 


A  +  r,,\  (i  -  r„\ 

\— — )sin  h  ~ 


(24) 


For  small  numerical  distances  (i.e.,  F(w)  -  1),  A  -  1,  and  B  -  sin  ^  -  A. 
Furthermore,  for  sin  ^ 1  »  |a|,  A  -  1  and  B  *  sin  ijij,  When  sin  is  compa¬ 
rable  to  or  less  than  A,  the  horizontal  distance  p  will  be  much  greater  than 
the  sum  of  the  transmitting  and  receiving  antenna  heights  (z  +  h)>,  In  the 
limit  as  ^  approaches  zero,  A  -  F(wQ)  and  B  -  -AF(wq),  where 


?K) 


1  -  i(ww  )1^2e”w°erfc(iw 


1/2 


(25) 


and 


wo  =  - 


vA 


(26) 


For  this  case  (i.e.,  p2  »  (z  +  h)2),  Wait*0,13  has  shown  that  F(w)  can 
be  replaced  by 


F (w)  -  [1  +  YqA(z  ♦  h)JF(w0) 

and  we  can  make  use  of  his  tabulated  results^  of  the  function  F(wQ). 


(27) 


Since  the  factor  A  (equation  (22))  is  different  from  unity  only  when  (1) 
the  angle  \|i,  is  very  small  and  (2)  the  Soramerfeld  attenuation  function  F(w) 
is  different  from  unity,  A  is  only  a  farfield  surface-wave  term.  Therefore, 
we  can  discard  all  derivatives  of  A  that  are  not  farfield  terms.  For  example, 
when  jn2|  >>  1  and  |a  sin  |  <<  1, 
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3  ( 

f  -YnR  \ 
Ae  P  1  j 

3p  \ 

(  R1  / 

-vrR,  ~Y«R, 

n  ^  ,  e  0  1  e01  3A 

-A(l  +  y  R  )cos  i|» - - —  + - x  — 

01  1  Rt  R,  3p 


-  -A(l  +  y0Rj)cos  t*/j 


-(!  +  YqRjAJcos  ifij 


kl  1 

,-YoRl 


Rl 


(28) 


>-^oRi 


R1 


Therefore,  from  equations  (25)  and  (28), 


$  .  j  .  .  r.  C_QS  t 


47rioje 


0 


'^0**0 


(1  +  YqRq)cos  ii<Q 


-  (1  +  yuR,)cos  ^  — -y-  + 


„  “Yo^i 

"YqRj  2  cos  ij^e  0  1 


(29) 


n2R2 


(1  +  YqRjA) 


which  is  identical  to  equation  (39j  of  Bannister.9 

rrom  equation  (11),  utilizing  the  identity  (Uj  -  uQ)(u1  +  uQ)  =  y2  -  Yq, 
we  have  v  ' 


ui  -  uo 


1 


Ylu0  f  ^0U1  uo(ui  +  V  U0(,J0  +  a2ui} 

Therefore,  equation  (2)  reduces  to 

„„-uoU+«)  r  ,-u0(z+h) 


(30) 


=  P  cos  4>  J_ 


-u0(z+n)  r  -uQ(z+h) 

— - -JQ(Ap)AdA  -  2A2  /  — H - 7? — 7  J0(Ap)AdA 

uo(ui  +  uo>  J0  uo(uo  +  A  ul) 

’■]  *  [ 


4iriue0  3p 
small  Sommerfeld  numerical 


(31) 


[distance  term 

Equation  (31)  is  equivalent  to 

_  p  cos  <j>  v  a 


term  to  account  for 
.  larger  numerical  distances 


]• 


n  =  r_.  z  x  Jt-(T  .  i  )  , 
2  4iriu>e0  3p  21  22 


(32) 


where 


'21 


/pn* 

■^o  \ui  +  uo / 


0(z+h)  ,  ,,  ,A  .. 

uo 


(33) 
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and 


Since 


and 


then, 


I 

22 


/-u0(z+h) 

— — - | - 

uo'un  +  *2uJ 


JQ  (Ao) AdX  , 


dC 


(z+h) 

J1 


dA  . 


Bannister7*8  has  shown  that 


(34) 


[55) 


(5) 


(36j 


-  -  ^[sin  f2e  y°Rz  -  (sin  -  Y0d)e  Y°  *]  ,  (37) 

which,  after  some  manipulation,  can  also  be  expressed  [utilizing  equation 
(A-8)]  as 

31  cos  *  e  Y°R?-  (1  +  Y0d)cos  Y°Rl 

— _  - £ - - - -  ,  (■>») 

3o  R2+d+z+h  Rj+z+h 

where  sin  ip j  =  (z  +  h)/Rj,  cos  ijjj  =  p/Rj,  sin  =  (d  +  z  +  h)/R2>  and 
cos  i|)2  =  p/R2- 

If  we  follow  the  same  procedure  as  in  the  derivation  of  V  •  n,  equation 
(34)  reduces  to 


1 


22 


-u  (z+h) 
e  u _ 

u0(u0  +  1 0 ^ 


J0(Xp)XdA 


=  2A2P  _  2P  _ 
Y0A  Yl 


(39) 
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and,  retaining  only  terms  in  1/R,  we  have 
31 


?? 


-  y„dP  cos  to  =  v  d  cos  to 
dp  0  i  o 


[1  -  F(w) ] 


,e 


-ynR 


0"1 


Therefore,  from  equations  (22),  (32),  (37),  and  (40), 

n7 


p  cos  j> 
4Triue0p 


,  ">0R2  ,  ‘^0R1 
sm  ip  e  0  z  -  sin  e  u  1 


+  Y0de  Y°  !(sin2  sJ/j  +  A  cos2  . 


(40) 


(41) 


Since  we  have  now  derived  expressions  for  the  HED  Hertz  vector  [equations 
(8),  (29),  and  (41)],  the  fields  in  air  can  be  obtained  from 


£  =  -Y2n  +  v(v  •  n) 

H  =  iu)CQ  (V  •  fl)  . 


(42) 


If  we  follow  the  same  procedure  as  outlined  above,  we  also  can  obtain 
suitable  expressions  for  the  HMD,  VED,  and  VMD  Hertz  vectors.  The  resulting 
HED,  HMD,  VED,  and  VMD  field-component  expressions  for  the  air-to-air  propa¬ 
gation  case  are  presented  in  tables  1  and  2.*  They  are  strictly  valid  for 
1 n2 |  >>  1.  However,  for  most  cases,  the  requirement  that  jn2j  >_  10  is  suf¬ 
ficient. 


These  formulas  reduce  to  the  author's  previously  derived  results  when 
either  (1)  the  Sommerfeld  numerical  distance  is  small, 7, 8  (2)  the  measurement 
distance  is  much  less  than  a  free-space  wavelength, 1 >2  or  (3)  the  measurement 
distance  is  much  greater  than  an  earth-skin  depth. 9  When  Oj  -*■  °°,  d  -*■  0, 

A  -+  0,  u  -*■  rl(  -*•  1 ,  R,  -*•  Rj,  F (w)  ->1,  A  -*■  1  and  the  formulas  reduce  to 
well-known  equations  for  propagation  over  a  perfectly  conducting  flat  earth. 
For  convenience,  these  equations  are  listed  in  tables  3  and  4. 

It  should  be  noted  that  the  VED  Ez  and  and  HMD  Ez  field-component 
expressions  presented  in  tables  1  and  2  were  first  derived  by  Norton. 14, 15 
Also,  all  VMD  components,  as  well  as  the  HED  and  HMD  Hz  components,  are  ident¬ 
ical  to  the  author's  previously  derived  results. ^>7, 8  por  the  sake  of  conven¬ 
ience,  we  have  tabulated  them  in  this  report.  In  these  tables  (tables  1 
through  4),  sin  ipQ  =  (z  -  h)/RQ ,  cos  i|>0  =  p/RQ,  sin  ^  =  (z  +  h)/R1 ,  cos  ^  = 

p/'Rj,  sin  <p2  =  ( d  +  z  +  h)/R2,  and  cos  ip2  =  p/R2. 


*A11  tables  have  been  placed  together  at  the  end  of  this  report. 
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SURFACE-TO-AIR  PROPAGATION 


The  HED,  HMD,  YED,  and  VMD  field-component  expressions  for  the  surface- 
to-air  propagation  case  (h  =  0,  z  0)  can  be  obtained  from  the  air-to-air 
propagation  equations  (tables  1  and  2)  simply  by  setting  h  =  0.  The  resulting 
equations  are  presented  in  tables  5  and  6.  In  these  tables,  R2  =  p2  +  z2, 

R?  =  p2  +  (d  +  z)2,  sin  p  -  z/R,  cos  ip  =  p/R,  sin  ip.  =  (d  +  z)/R.,  and  cos  i. 

1  ill 

=  P/R^.  These  formulas  are  strictly  valid  for  |n2|  >>  1.  However,  for  most 
cases,  the  requirement  that  |n2j  >_  10  is  sufficient. 

Image-theory  expressions  for  the  subsurface-to-air  propagation  case  can 
be  obtained  from  the  surface-to-air  propagation  equations  (tables  5  and  6) 
simply  by  multiplying  each  expression  by  exp(y,h).  fAll  VED  components  must 
also  be  multiplied  by  1/n2  to  satisfy  the  boundary  conditions.)  The  resulting 
formulas  will  be  valid  for  j  n2i  »  1  and  R2  »  |  h  | 2  .• 


AIR-TO-SURFACE  PROPAGATION 

The  HED,  HMD,  VED,  and  VMD  field-component  expressions  for  the  air-to- 
surface  propagation  case  (h  ^  0,  z  =  0)  can  be  obtained  from  the  air-to-air 
propagation  equations  (tables  1  and  2)  simply  by  setting  z  =  0..  The  resulting 
equations  are  presented  in  tables  7  and  8.,  In  these  tables  D2  =  p2  +  h2, 

D?  =  p2  +  (d  +  h)2,  sin  !|p  =  h/D,  cos  it*  =  d/D,  sin  =  (d  +  h)/Di,  and  cos 
=  p/D^ .  These  formulas  are  strictly  valid  for  |n2)  >>  1.  However,  for  most 
cases,  the  requirement  that  J n2  j  ^  10  is  sufficient. 

Image-theory  expressions  for  the  air-to-subsurface  propagation  case  can 
be  obtained  from  the  ?i r-to-surface  propagation  equations  (tables  7  and  8) 
simply  by  multiplying  each  expression  by  expfyjZ).  (All  Ez  components  must 
also  be  multiplied  by  1/n2  to  satisfy  the  boundary  conditions.)  The  resulting 
equations  will  be  valid  for  |n2  |  »  1  and  D2  >>  |z|2. 


SURFACE-TO-SURFACE  PROPAGATION 


The  (simple  form)  HED,  HMD,  VED,  and  VMD  field-component  expressions  for 
the  surface-to-surface  propagation  case  can  be  obtained  from  the  air-to-air 
propagation  equations  (tables  1  and  2)  simply  by  setting  both  z  and  h  equal  to 
zero.  The  resulting  equations  are  listed  in  table  9.  In  these  tables,  p?  = 

P2  +  d2,  sin  =  d/p^,  3nd  cos  =  p/p^.  These  formulas  are  strictly  valid 
for  |n2(  »  1.  However,  for  most  cases,  the  requirement,  that  jn2|  ^10  is 
sufficient. 

It  should  be  noted  that  the  VMD  Hp  (and,  by  reciprocity,  the  HMD  Hz) 
image-theory  expressions  for  the  surface-to-surface  propagation  case  are  only 
valid  when  p  »  6,  where  6  is  the  earth-skin  depth.  When  ‘.his  condition  is 
not  satisfied,  other  formulas  should  be  utilized  for  these  cwo  field  components 
(for  example,  see  table  9  of  Bannisterl6) . 
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Image-theory  expressions  for  the  subsurface-to-subsurface  propagation 
case  can  be  obtained  from  tiie  surface-to-surface  propagation  equations  (table 
9)  simply  by  multiplying  each  expression  by  exp[Yj(z  +  h)].  (The  "ED  Ep  and 
and  the  HED  and  HMD  E„  components  must  be  multiplied  by  l/n°,  while  the 
VED  E„  component  must  be  multiplied  by  1/n4,  to  satisfy  the  boundary  condi¬ 
tions.)  The  resulting  equations  will  be  valid  for  i  n‘?  i  >:>  1  and  p2  >> 

|z  +  h]2.  It  should  be  noted,  however,  that  more  accurate  formulas  are  avail 
able  for  the  subsurface-to-subsurface  propagation  case  when  ! n2 I  >>  1 . 16 , 17 


CONCLUSIONS 


In  this  report,  we  have  extended  the  use  of  finitely  conducting  earth- 
image  theory  techniques  to  any  range  and  have  derived  formulas  for  the  elec¬ 
tric  and  magnetic  fields  produced  by  the  four  elementary  dipole  antennas  for 
the  air-to-air,  surface-to-air,  air-to-surface,  and  surface-to-surface  propa¬ 
gation  cases.  The  only  restriction  on  the  use  of  these  formulas  is  that 
jn2|  >_  10.  They  are  valid  at  any  frequency  and  at  any  range  for  the  flat- 
earth  case.  These  formulas  reduce  to  the  author's  previously  derived  results 
when  either  (1)  the  Sommerfeld  numerical  distance  is  small,  (2)  the  measure¬ 
ment  distance  is  much  less  than  a  free-space  wavelength,  or  (5)  the  measure¬ 
ment  distance  is  much  greater  than  an  earth-skin  depth.- 


The  results  presented  in  this  report  can  be  extended  to  a  multilayered 
earth  simply  by  letting  d  =  (2/yj)Q,  where  0  is  the  familiar  plane-wave  cor¬ 
rection  factor  employed  to  account  fc.r  the  presence  of  stratification  in  the 
earth. 13  For  a  homogeneous  ground,  the  argument  of  the  numerical  distance  wQ 
is  always  between  0  and  -90,  resulting  in  the  transverse  magnetic  (TM)  sur¬ 
face-wave  fields  varying  as  1/p2  as  p  -*■  For  a  stratified  ground,  the 
argument  of  wQ  can  be  positive,  resulting  in  the  TM  surface-vave  fields 
varying  as  l/^p\13 


It  should  be  noted  that  the  two  media  can  be  inverted  and  the  air 
replaced  by  the  earth's  crust  (of  conductivity  a2  and  dielectric  constant  e 
The  same  equations  (tables  i  through  9)  can  be  utilized  as  long  as  |n2|  = 


Y2/y2  I  >_  10  simply  by  replacing  iu>e0  by  a?  +  iwe-,. 


These  formulas  are  intended  to  supplement  the  author's  recently 
derived^ >16  subsurface-to-subsurface  and  a;r-to-air  propagation  formulas,  in 
terms  of  computer  time,  the  formulas  presented  in  these  three  reports  can  be 
evaluated  in  fractions  of  a  minute  compared  with  hours  fer  the  complete  numer¬ 
ical  evaluation  of  the  exact  Sommerfeld  integrals. 
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,‘Y0*0 


(3  ♦  3y0ro  *  YoRo>sin  *0  cos  *0  — p- 


4viueQ 
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VED 
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S~Y0R1 
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e'YoRi 
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|c-Field  Air-to-Air  Propagation  Formulas  (|n2j  :>  10) 

t 


Table  1.  (Cont'd)  Electric-F 


Dipole 

Type 


HED 


p  cos  ♦ 


4«iuc 


0 


(3  cos2  ♦q  -!)(!♦  y0R0)  -  YqRo  si"2  *0 


-Y„R 


0  0 


r3 

K0 


(3  cos2  -!)(!♦  y0R!)  -  r„YgRf  sin2  *  ' 


*YnR 


0  1 


■5f 


-y.R, 
2  e  3  1 


2Ri 


K3  ~s2  tx  -  l)Ci  ♦  *,«,)]  -  Tz 


iAe’VVV 

R, 


♦  n2y2R2 


sin  ^  ^ — ~^AF(w)j  J 


P  sin  » 


4ai»c. 


-  (1  ♦  Y 


-Y„R 

2e  ® 


ft"Rj 


>MD 


iwpQa  cos  4 


4* 


.-*0R0 


O  ♦  YQRj)sii>  4q  3T 
*o 


-t0ri  /l  -  r  \  -y0ri 

^  *  rnYoRi)sin  -  yF(w)y^R2  jqf 


-y„R 


0  2 


(1  ♦  Y0d)e 


‘W 


R2(R2  ♦  d  ♦  z  *  h)  Rj(Rj  ♦  z  ♦  h) 


1MW0«  Sll 


4* 


♦  (1  ♦  Y, 


(R2  * 
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Table  2.  Magnetic-Field 


Dipole 

Type 


VED 


JL 

4w 


C1  +  Y0R0)cos 


♦  (1  -  rj,)F(w)cos 


VMD 


■  f  1  o  e"Y°R° 

—  1(3  ♦  3y0R0  ♦  YgRg)sin  *0  cos 

,  ,  e"Y0R2 

-  (3  ♦  3Y(JR2  ♦  Y^sin  ♦  cos  *2  — — 

*2  J 


4* 


-y„R 


(1  ♦  Y0R2)sin  *2  r2 


02 
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e‘Y0R0 

-  f1  ♦  Y0Ro)sin  ♦„  — -j— 

*0 


e'Y0R2 


— Y  ftRi 

(1  ♦  YndA)e  0  1 


R2(R2  ♦  d  ♦  z  ♦  h)  Rj(Rj  ♦  2  ♦  h) 


P.cos  f 


4* 


(1  ♦  Tfl 


’  t1  ♦  r.»Vl)sin 

-Y.R, 


'0  2 
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4* 


■Y0Ro 
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Ko 
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e'Y°R2 
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tMD 


e'T0Rl 
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”Y0R1  ] 

♦  1(2  ♦  2y0R,A)  -  sin2  ♦1(3  ♦  Sr^  ♦  Y^R2)]^—  J 


■  cos  4 


4* 


(1  ♦  Y0! 


♦  tr„  ♦  (i  -  r„)F 
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-  ^IU  ♦  YqR0  ♦  Yoro)  -  sin2  M3  4  3*0Ro  4  *oRoJl 
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Table  4.  Magnetic-Fi 
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Type 
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-vnR 


0  1 


R3 


HED 


_  p  sin  ♦ 


4w 


-yn*> 


(1  +  YgRpisin  <>0  J^2 


0  0 


,"Y0R1 


-  (1  ♦  Y0Rj)sin  — p- 


H4D 


■  sin  ♦ 


4w 


K2  ♦  2Y0 V  ’  sin2  *0(3  *  3y0R0  + 


♦  1(2  ♦  2y0R1)  -  sin2  ♦jO  ♦  SYpRj  ♦  Y§Rf )  ] 


c 

I 


I 
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Table  5.  Electric-Field  Surface-to-Air  Propagai 


-  *^--y[(l  -  3  sin2  \|i)(l  +  y„R)  +  Yqr2a  cos2  +] 
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Tab?'  6.  Magnetic-Field  Surface-to-Air  Propagation  F( 


Dipole 

Type 


_  -*oR 

—  (3  ♦  3v„R  ♦  Y^R2)sin  ^  cos  ^  — 5 — 
4ir  v  T0  *0  1  v  R3 


(3  ♦  3y0R.  +  Y2R?)sin  *.  cos  *. 


p  sin  ♦  „  ..  .  e  Y°Ri 

c 1  (1  ♦  Y«R  )sm  <>.  — •$ — 
4w  v  0  i  i  R. 

i 


-  U  ♦  Y0R)sin  ?  -p 


(1  ♦  Y0dA)e 
R(R  ♦  z) 


.‘YcRi 


R.  (R.  ♦  d  ♦  z) 

l  i 


[1  ♦  Y„d)e 


R(R  ♦  z) 


-  ~;±  ~T~[2  ♦  Y0R(1  ♦  A)  -  sin2  *(3  ♦  3y0R  ♦  y2R2)] 

-y0Ri 

♦  [(2  ♦  2y„R.  ♦  Y2R?)  -  Sir.2  *.(3  ♦  3y„R.  ♦  Y»R?)j  ■  r— 

0  1  0  1  lv  0  1  '0  1  J  R3 

-y0r 

-  1(2  ♦  2y0R  +  y|R2)  -  sin2  *(3  ♦  3yqR  ♦  YpR2)]"^- 


; 

»to-Air  Propagation  Formulas  (|n2|  >_  10),  [R2  =  p2  +  z2,  R?  =  p2  (d  ♦  z;2] 
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Dipole 

Type 


y,p  cos  ie  u 

- ="  YnDA  ♦  51— (3  ♦  3y«D  ♦  Y?D2) 

2ir(Oi  ♦  1(i)Ej)D2  0  YjD  0  0 

D3e^0P  (1  ♦  Y0d)e~Y°D  e'Y0Pi  _  D 

d  D(D  ♦  h)  Di(Di  +  d  +  h)  *  Y° 


rV  +  •  _  in 3  1^  COs2  *  '  +  Y0D) 

2*{o1  ♦  lwe^D3  I  0 

»  m 

-  1  '  5^  e  Y°^Dl  ^  *  nYgR2(sin  *  -  B) 


(1  +  Y^OJcos 


2it(0j  ♦  iuCjJD3 


d2  “i 


Y,«  cos  #e 


iw(jQ»  sin  $ 


-  (1  ♦  Y0»)sir. 


(1  ♦  Y()dA)e 


D(D  ♦  h) 


face  Propagation  Formulas  (,;n*.  >  101,  IP*  = 


h*.  P*  =  c*  ♦  *  hv 
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Table  8.  Magnetic-Field  Air-to-Surface  Propagi 


Dipole 

Type 


-  ^  (3  ♦  3yqD  ♦  Y|jD2)sin  1>  cos  *  ^ 


♦  (3  ♦  StpDj  ♦  YgD?)sin  cos 


LJ±±  (1  ♦  Y0D)sin  ♦ 


e‘Y0Di 

*  (1  ♦  -oD^sin 


(1  ♦  Y„dA)e 


D(D  ♦  h) 


Di(Di  ♦  d  ♦  h) 


-YqD 

"  S™~*  -"a-  -[2  ♦  yqD(1  ♦  A)  -  sin2  *(3  ♦  3yqD  ♦  y^D2)] 


♦  [(2  ♦  2y  D.  ♦  y^D?)  -  sin2  *.  (3  ♦  3y  D.  ♦  y2D?)1— 5- 
0  1  0  1  1  0  1  0  l  D? 


-r0Dl 

[(2  ♦  2y0D  ♦  y^D2)  -  sin2  t(3  ♦  3yqD  ♦  Y§D2)  3^3— | 


Cl  ♦  1 


a  cos 

4w 
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e  Propagation  Formulas  C | n2 |  >10),  [D2  =  p2  ♦  h2,  D2  =  p2  ♦  (d  +  h)2] 


[(1  ♦  Y0D  ♦  YqD2)  -  sin2  \l>( 3  +  3y„D  +  Y„D2)] 


0  'o  d3 


„-Y0Di 

-  K1  +  Y0Di  V  y^D2)  -  sin2  *.(3  ♦  3yQ[)i  -  YgD?)j . ^  - 


D(D  ♦  h)  Di(D<  +  d  ♦  h) 


d  ♦  Y0D  ♦  2AYqD2)~-|  ♦  (l  ♦  Yod.)! 


sin  6  e  u 

-^-Z-  (1  ♦  y0D)cos  tl>  -  (1  +  Y0Di)cos  *i 


_Y(jDi  "j  -  f  p  +  3y0D  ♦  YpD2)sin  *  cos  f 


'Y0Di 


(3  ♦  3y0D.  ♦  YjD2)sin  *.  cos  *.  — 
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Appendix 


USEFUL  APPROXIMATIONS  WHEN  R  »  |dj 


—  -  1 

R, 


d2  I  2Rj  sin  *,  \ 

iSfV* - 3 - 3sln  *-) 


G:j 


i  sin 


-  4  sin2  ip 


/*  lY  .  3d2/,  2R,si">'l  r 

\rj  ' 1  -  HfV1 - 1 - 51  * 


.) 
.) 

/*lY  ,  3d2/  2R,  si"  *,  ,  .  ,  .  \ 

I  — I  -  1  -  -  t-II  +  - - - 7  sin2  i I 

\  R2/  2Rf\  d  7 


6  sin2  iji 


and 


(sin  d  cos2  \p  \ 

— - — j 

For  1 7  x  C  —  +  h) |  »  1  (i.e.,  |2(R]/d)sin  »  1) , 

/  R  V  xd  sin  ^ 

\rT/  " 1  ‘ 


72R2d2  sin2 


2R2 


and 


/WV  , 


1  -  Y 


(d  sin  <i>  \ 

—)■ 


Some  other  useful  approximations  are 

-ynR,  -ynR, 

cos  4 i^e  u  2  (1  +  YQd)cos  j e  0  1 


(A—  1) 


(A-2) 


(A-3) 


(A-4) 


(A-S) 


(A-6) 


(A-7) 


(A-8) 


R2  +  d  +  z  +  h 


+  z  +  h 


e’Y°Rl 

'  -  d  cos  *l(l  +  Y^)— -r~ 


(A-9) 


A-i 


^  .ariii  ,i  - 
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_Y  nR 


0  2 


(1  +  -Y0d)e 


"Y0R1 


R  (R  +  d  +  z  +  hj  RjCRj  +  z  +  h) 


,“Y0R1 


'  d(1  +  W*  13 


and 


=.  Y  0R2 


(1  +  YQdA)e 


'0  '1 


1 

Y0R' 


R2(R2  +  d  +  z  +  h)  R,(R,  +  2  +  h)  J 
R 


2R2 

7^ 


! _ L  e  Yo  (R2  ri) 

Ro 


'  (1  +  y,Rj  sin  i|0(l  +  Y0Rj) 
-  sin2-  p 1  (3  +  3y 0 R j  +  YqR^) 


,R„  -VnR, 

-  -  (1  +  YgR. J  sin  ^ 


e'YCK2  e*Y0Kl 

(1  +  >0R2)sin  *0  — p- 


R 


de 


-Y  aR 


0*  1 


3  1(1  +  Y0 Rj)  -  sin2  ^(3  ♦  3YqR1  +  YqRjJ 

K1 


e^*2 . e'YoRi 


(1  +  y0R2)cos  v2  — p—  -  (1  +  V^05  p: 


Rf 


d  sin  iji  cos  1  e 


Y0R1 


R? 


-(3  ♦  3y0R1  -  Tr|R{)  • 


(A- 10) 


+  Y0RlA-)  rT—  i-a'11) 


(A-l?) 


CA- 13) 


(A-14) 


where 


R1 = 

:  D2 

+  (2  + 

h);'. 

R2  = 

n 

-i 

■  p- 

+  (d  + 

z  +  h) 2 

sin 

*1 

=  (2  * 

h)/Rlt 

cos 

Vl 

=  p/Rj, 

sin 

V2 

=  (d  + 

z  +  h)/R 

cos 

*2 

=  p/R2. 
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